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I .  INTRODUCTION 


During  che  past  year,  Che  research'  carried  out  under  the  grant  was 
directed  to  two  areas:  (1)  the  development  of  microstructure  in  thermo¬ 
setting  resins,  most  notably  epoxy  resins;  and  (2)  the  exploration  of  the 
chemistry  and  kinetic  processes  involved  in  the  conversion  of  organo- 
metallic  materials  to  ceramics.  Progress  in  each  of  these  areas  will  be 
summarized  in  the  sections  which  follow. 

The  student  working  on  the  development  of  microstructure  in  thermo¬ 
setting  resins  will  complete  his  doctoral  thesis  in  September.  His  place 
will  be  taken  by  a  second  graduate  student  investigating  the  conversion 
of  organo-metallic  precursors  to  ceramics.  As  indicated  in  the  appended 
proposal,  the  latter  area  will  represent  the  exclusive  focus  of  our  work 
for  the  coming  year. 


-Y 

II.  MICROSTRUCTURE  OF  EPOXY  RESINS 

t- 

\  i  . The  present'work  represented  a  continuation  of  our  efforts  directed  to 
understanding  the  microstructure  of  epoxy  resin  cured  with  amine  hardeners. 

The  focus  has  been  on  characterizing  a  series  of  networks  formed  by  changing 
the  formulations  and  the  processing  conditions  of  the  system.  The  mor¬ 
phological  analysis  has  been  supplemented  by  a  statistical  model  for  the 
network  formation;  and  the  results  have  been  used  to  understand  and  control 

the  mechanical  properties  of  these  materials. 

a 

The  table  on  the-Yollowing  pagelshows  the  Epon  828  epoxy/ triethylene 
tetramine  (TETA)  formulations  which  were  explored,  together  with  the  curing 
conditions  employed. 

Figure  1  shows  the  glass  transition  temperatures  (T  )  of  these  specimens. 

8 

Table  2  shows  the  molecular  weight  between  the  crosslinks  (M  )  for  some  of 

c 

the  formulations.  was  evaluated  from  the  rubbery  modulus  of  the  resins 
after  applying  rubber  elasticity  theory.  The  epoxy-rich  formulations 
have  higher  M  ,  which  decreases  towards  the  stoichiometric  formulations, 
and  then  again  increases  with  increasing  curing  agent  concentration.  M 

c 

decreased  after  the  post  curing  treatment.  A  relationship  has  been  found 

between  T  and  M  : 

g  c 

T  =  51.7  +  23800/M  (1) 

g  c 

The  modulus  at  room  temperature  of  these  resins  is  shown  in  Figure  2. 

The  modulus  decreases  with  increasing  hardener  content,  especially  around 
stoichiometry.  From  the  published  literature,  epoxies  have  been  assigned 
a  two  phase  structure  (nodular)  from  an  analysis  of  their  fracture  surfaces 
and  dynamic  mechanical  properties.  The  dispersed  phase  accounts  for  more 


3. 


TABLE  1 

EPON/TETA  Formulations 


Condition 

Concentration  of 
Curing  Agent ,PHR 

Epoxide/Amine 
Hydrogen  (E/A) 

Curing  Condition 
hour  ,  C 

A 

6.5 

2:1 

7,120 

B 

6.5 

2:1 

4  RT/7,120 

C 

8 

1.63:1 

4  RT/7,120 

E 

10 

1.3:1 

7,120 

F 

10 

1.3:1 

4  RT/3,83 

G 

10 

1.3:1 

4  RT/8,83 

H 

12 

1.06:1 

1  RT/8,83 

11 

13 

1:1 

7,110 

12 

13 

1:1 

7,110 

K1 

13 

1:1 

8,83 

K2 

13 

1:1 

8,83 

LI 

13 

1:1 

1  RT/8,83 

L2 

13 

1:1 

1  RT/8,83 

M 

17 

1:1.3 

12,80 

N 

17 

1:1.3 

3  RT/9,80 

Q 

26 

1:2 

12,80 

R 

26 

1:2 

3  RT/9,80 

*RT  is  Room  Temperature. 


A  part  of  each  of  the  samples  was  also  subjected  to  Postcure  at  150  C 
for  two  hours.  This  will  be  indicated  in  the  text  by  a  superscript  P. 


EPON/TETA 
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7. 


than  50  volume  %  with  an  average  nodule  size  of  200-300  A.  The  nodules  are 
assumed  to  be  rich  in  amine  and  thus  have  a  higher  crosslink  density  than 
the  matrix.  This  increases  their  hardness  and  makes  them  undeformable. 

Figure  2  clearly  shows,  however,  that  amine-rich  formulations  have 
lower  modulus  than  epoxv-rich  or  stoichiometric  compositions.  Hence  the 
arguments  based  on  the  occurrence  and  characteristics  of  nodules  are  incon¬ 
sistent  with  the  data.  It  has  also  been  shown  by  SAXS  studies  that  if 
nodules  did  exist  in  the  sizes  and  volume  fractions  claimed  above,  then 
the  density  difference  between  the  matrix  and  nodules  would  be  less  than 
0.1%. 


Figure  3  shows  a  crossplot  of  the  T  vs  modulus  data.  At  both  the 

8 

ends  of  the  curve,  T  drops  without  a  change  in  modulus.  On  the  left  hand 

S 

side  of  the  figure,  the  drop  occurs  for  materials  which  are  amine  rich 
(conditions  Q,  R,  QP,  RP)  and  at  the  other  end  for  the  epoxy  rich  materials 
(A,  B,  AP,  Bp).  This  drop  cannot  be  taken  as  the  formation  of  the  second 
phase,  because  if  this  happened  then  the  modulus  would  change  at  constant 


This  phenomenon  can  be  explained  by  a  plasticization  effect.  In  the 
epoxy-rich  region,  unreacted  epoxide  plasticizes  the  matrix;  and  in  the 
amine-rich  region,  the  high  concentration  of  the  flexible  curing  agent 
gives  the  same  effect.  STEM  observations  (under  Z  contrast)  on  microtomed 
and  stained  sections  of  these  specimens  show  a  heterogeneous  nature  (Figure 
A).  These  were  observed  only  for  amine  rich  formulations,  or  for 
stoichiometric  compositions  which  were  precured  at  room  temperature.  DSC 


and  DMA  do  not  show  the  presence  of  two  T  's  or  shoulders  in  Tan  6  around  T  . 

8  8 
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Figure  3  -Onset  of  Glass  Transition  Temperature  from  DSC  vs  Flexure  Modulus. 


9. 


Figure  A  -  Stained  STEM  micrographs  for  stoichiometric 

(condition  L2P)  formulation.  Room  temperature 
cure  was  given. 
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It  has  previously  been  observed  in  an  aromatic  ami ne/DGEBA  epoxv  svstem 
that  the  reactivity  ratio  of  the  primary  and  secondary  amine  (K^/K,,)  is 
around  5.  If  this  were  the  case  in  the  present  system,  then  in  an  amine 
rich  system  the  primary  H  will  react  with  the  epoxies  tending  to  form  long 
linear  chains,  which  could  then  precipitate  out  of  the  epoxy-rich  matrix 
creating  a  two-phase  system. 

It  was  decided  therefore  to  model  the  growth  of  a  system  with  a  varying 
K.^/K-2  and  stoichiometry  to  obtain  insight  into  the  types  of  species  formed 
at  various  stages  of  reaction.  Figure  5  show  the  fractions  of  primary  (Ql) 
amine  reacted  for  different  stoichiometries  and  reactivity  ratio  r(r=K^/K0) 
as  a  function  of  epoxide  consumption.  This  figure  shows  the  variations 
only  up  to  the  onset  of  gelation.  It  should  be  noted  from  the  figure  that 
the  epoxide  consumed  at  the  gel  point  does  not  vary  greatly  with  changing 
reactivity  ratio. 

Table  3  shows  the  various  species  that  can  be  generated  during  the 
curve  of  epoxies;  and  Tables  4  to  9  show  the  concentrations  of  the  different 
components  at  various  stages  of  the  reaction.  It  is  seen  that  at  the  gel 
point,  in  the  epoxy  rich  mixture,  most  of  the  contribution  to  the  total 
extent  of  reaction  comes  from  highly  branched  units  (N2PIS,  N2P2S,  N2p3S); 
whereas  in  amine-rich  mixtures  (A/E=2),  it  comes  from  the  less  branched 
units  (N1P0S,  N1P1S,  N0P1S,  N2P0S).  The  only  branched  units  contributing 
substantially  to  the  latter  are  N2P1S  or  N1P2S. 

With  increasing  reactivity  ratio  (K^/K^),  one  observes  that  in  amine 
rich  mixtures  the  amount  of  unreacted  TETA  is  reduced  and  the  contribution 
from  N1P2S  decreases  at  the  expense  of  N2P1S  (although  both  have  nree  of 
their  functionalities  reacted). 


7.  Epoxide  Reacted 

Fraction  of  primary  amine  reacted  with  increasing  consumption 
results  are  plotted  either  till  gel  point  or  till  the  primary 
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TABLE  3 

Types  of  Different  Species*  Generated  for  TETA. 


NOPQS 

pssp 

N1P0S 

psss 

N2P0S  , 

SSSS 

N0P1S 

PSTP  ‘yJ — 

N0P2S 

PTTP 

N1P1S 

PSTSy  PTSS,  TSSP 

N1P2S 

PTTS t  TSTP,  TTSP 

N1P3S 

TTTP 

N2P1S 

STSSy  SSST 

N2P2S 

N2P3S 

STTS,  TSSTy  TSTS,  TTSS 

TTTS,  TSTT 

N2P4S 

tttt 

*  Specie  N0P1S  indicates  zero  primary  and  one  secondary  have 
reacted.  This  can  also  be  written  as  PSTP,  i.e.,  the  first 
nitrogen  is  primary,  the  second  is  secondary,  the  third  is 
tertiary  and  the  last  one  is  also  primary.  The  schematics 
show  TETA  molecules,  a  solid  end  indicates  -that  it  is  reacted. 


TABLE  4 

%  of  Different  Species  at  Increasing  Epoxide  Conversion  (EPON  /  TETA) 


TABLE 


TABLE  8 

X  of  Different  Species  at  Increasing  Epoxide  Conversion  (EPON  /  TETA) 
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TABLE  9 

X  of  Different  Species  it  Increasing  Epoxide  Conversion  (EPON  /  TETA) 
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To  examine  whether  such  results  for  large  r  (>2)  in  amine-rich  systems 
would  lead  to  the  formation  of  long  linear  polymer  chains  (which  would  then 
precipitate  to  give  a  second  phase),  further  calculations  were  carried  out 
^see  appendix  for  details  of  calculations),  with  the  results  shown  in 
Tables  10  and  11.  These  tables  lists  the  %  probability  for  the  formation 
of  diads  and  triads.  The  type  of  triad  listed  is  of  the  type  which  will 
lead  to  longer  polymeric  species.  The  probability  of  the  formation  of  this 
triad  increases  appreciably  at  r  =  5  in  amine  rich  formulations.  In  the 
limit  that  r  =  2,  the  reactivity  differences  between  primary  and  secondary 
amines  will  not  cause  an  appreciable  change  in  the  nature  of  the  species. 
However,  it  may  help  explain  the  inhomogeneity  of  amine  rich  mixtures. 

To  consider  this,  it  is  useful  to  evaluate  the  difference  in  the 
solubility  parameters  required  for  precipitation.  It  has  been  argued  that 
precipitation  occurs  when  A  is  greater  than  5 


for  polymeric  species,  and  greater  than  10  for  monomeric  species.  6^  and 
6^  are  the  contributions  to  the  solubility  parameters  due  to  dispersive 
polar  forces  and  hydrogen  bonding  respectively.  Subscripts  1  &  2  refer 
to  the  precipitating  species  and  the  matrix  component.  Table  12  lists 
values  of  these  parameters  for  different  species.  It  can  be  argued  that 
the  polymeric  species  shown  in  Table  13  has  a  A  of  9.28  (as  compared 
to  DGEBA)  which  is  more  than  5  hence  it  would  precipitate.  The  reference 
point  has  been  taken  as  DGEBA  because  only  43%  of  epoxide  had  reacted  at 
gelation;  therefore  an  assumption  was  made  that  a  sea  of  unreacted  epoxy 


remained. 
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TABLE  10 

Percentage  Probability  for  the  Formation  of  Triads  as  the  Function  of 
Stoichiometry  and  K1/K2  at  different  Stages  of  P.eaction 

%  Primary  Amine  Reacted 


Ki/K2 

15 

A/E  = 

30 

2 

Gel 

30% 

A/E  =  1 

60% 

Gel 

30% 

A/E  =  0. 

60% 

,5 

Gel 

1 

8.99 

15.64 

8.47 

7.81 

2.71 

1.08 

3.9 

1.36 

0.022 

2 

9.09 

19.93 

10.4 

9.97 

6.95 

0.94 

4.99 

3.46 

0 

5 

8.9 

22.3 

7.12 

11.14 

11.46 

0.12 

5.57 

5.74 

- 

7 


TABLF.  11 

* 

Percentage  Probability  for  the  Formation  of  Diads  as  the  Function  of 
Stoichiometry  and  K^/K.2  at  Different  Stages  of  Reaction 

%  Primary  Amine  Reacted 


A/E  «  2 

A/E  =  1 

A/E  = 

0.5 

Kx/K2 

15 

30 

Gel 

30 

60 

Gel 

30 

60 

Gel 

1 

26.3 

27.0 

10.4 

33.6 

14.1 

7.56 

36.9 

18.2 

.1.7 

2 

24.9 

30.2 

9.4 

35.7 

22.5 

5.75 

38.4 

27.2 

0.6 

5 


23.8 


32.2  6.64 


37.0  30.8  1.9 


39.3  35.6 


TABLE  12 

Solubility  Parameters  of  Different  Species  (( J/cm3 ) 1 /2 


Type  of  Species 

*d 

4p 

6v 

Sh 

Triad 

16.44 

1.743 

16.53 

11.61 

Primary  Diad 

17.3 

2.04 

17.42 

10.06 

Secondary  Diad 

17.25 

2.60 

17.44 

10.81 

DGEBA 

17.21 

0.64 

17.22 

1.72 

*  ov  *  '/|5d2  +  6p2 


Subscripts  d,  p,  and  h  refer  respectively  to  dispersive,  polar  and 
hydrogen  bonding  contributions  to  the  solubility  parameter. 


TABLE  13 


Difference 

JL 

in  Solubility 

(A)  for  Different 

Species 

Species  1 

Species  2 

DGEBA 

Diad  Primary 

Diad  Secondary 

DGEBA 

0 

8.34 

9.09 

Diad  Secondary 

9.09 

0.75 

0 

Biad  Primary 

8.34 

0 

0.75 

Triad 

9.91 

1.79 

1.21 

*  A  =  /  (<5vl  -  <5v2)2  +  ^hl  “  &h2^2 


From  a  probabilistic  calculation,  the  completely  undreacted  molecules 

2 

of  epoxide  will  be  only  (1-43)“ — i.e.,  32%.  At  gelation  (gelation  occurs 
at  a  different  value  in  amine-rich  system  here  due  to  a  different  stoichio¬ 
metry),  completely  uncreacted  epoxy  is  present  in  only  a  small  quantity 
which  cannot  form  the  matrix  (Table  14).  The  matrix  consists  of  diads  and 
other  species  which  have  solubility  parameters  (Tables  12,  13)  intermediate 
between  those  of  unreacted  epoxy  and  the  triads  (or  the  higher  polymeric 
species  that  are  formed).  Considering  one  step  before  gelation,  where  the 
phase  separation  can  take  place  without  being  hindered  sterically,  one  can 
consider  that  the  diads  (wnich  are  still  soluble  in  epoxy)  act  as  buffers  to 
keep  the  triads  in  solution.  The  difference  in  solubility  parameter  between 
the  epoxy  and  the  triads  can  be  reduced;  and  it  is  concluded  that  phase 
separation  giving  rise  to  a  large  fraction  (50  volume  %  or  more)  of  hetero¬ 
geneities  cannot  be  explained  in  this  way. 


TABLE  14 


Various  Kinds  of  Epoxide  Present  (%)  at  Different  Stages  of 
Reaction  as  a  Function  of  Stoichiometry  and  K^/K2 


K1/K2  Type  of 
Species* 


« - — • 


X  of  Primary  Amine  Reacted 


A/E  = 

2 

A/E  = 

1 

A/E  « 

0.5 

Gel 

30 

15 

___ 

60 

30 

Gel 

60 

30 

13.5 

•  5.3 

71.2 

30. A 

A0 

70 

A6.8 

66.6 

84 

Af  .5 

AA 

26.3 

A9.5 

A6.5 

27.3 

A3. 2 

30 

15 

AO 

10.7 

2.43 

20.1 

13.5 

2.7 

10 

3A 

1.0 

12.0 

53.9 

75.9 

29. A 

50 

75.2 

A6 

72.9 

87 

A5.3 

39 

22. A 

A9.6 

Al.A 

23.1 

A3. 6 

25 

12. A 

A2.7 

7.1  ' 

1.65 

21 

8.6 

1.77 

10.3 

2.1 

0.4A 

9.62 

59.7 

78.0 

27.7 

57.9 

78.6 

- 

77.5 

89 

A2.8 

35.1 

19.8 

A9.9 

36. A 

20.1 

- 

21.1 

10.7 

A7.6 

52 

1.25 

22. A 

5.73 

1.29 

— 

1A.3 

0.3 

*  A  solid  end  indicates  reacted. 


one  end  reacted. 


shows  an  epoxide  molecule  with 
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'll!.  CERAMICS  FROM  ORGANO-METALLIC  PRECURSORS 

Our  investigations  ot  the  conversion  of  organo-metallic  precursors  " 

to  ceramics  have  directed  attention  to  five  principal  areas.  These  were: 

1.  Investigation  of  hydrolysis  and  condensation  kinetics  of  tetra- 
ethylorthosilicate  (TEOS) ,  an  important  precursor  for  many  glasses  and 
ceramics ; 

2.  Exploration  of  the  competition  between  crystallization  and  viscous 
sintering,  to  determine  the  conditions  under  which  dense  glassy  bodies  or 
coatings  can  be  obtained; 

3.  Investigation  of  the  hydrolysis  and  condensation  of  tantalum 
ethoxide  to  form  coatings  of  tantalum  oxide; 

4.  Initial  exploration  of  organic  modified  oxides  as  novel  coatings 
and  bulk  materials;  and 

5.  Exploration  of  the  conversion  of  aluminum  alkoxides,  most  notably 
aluminum  sec  butoxide,  to  aluminum  oxide. 

Progress  in  the  first  two  of  these  areas  will  be  described  in  the  present 
report;  progress  in  the  last  three  will  be  deferred  to  a  later  date. 

/N 

A.  Hydrolysis  and  Condensation  Kinetics  of  Tetraethylorthosilicate  (TEOS) 

The  sol-gel  process  for  preparing  glasses  and  crystalline 
ceramics  has  many  potential  advantages  over  conventional  melt  or  powder 
processing  techniques.  Among  these,  four  have  specific  relevance  to  the 
present  work.  These  are:  (1)  the  high  level  of  chemical  homogeneity  which 
is  attainable;  (2)  the  possibility  of  relatively  low- temperature  processing; 

(3)  the  ability  to  tailor  composition,  structure  and  properties;  and  (4) 
the  ease  and  cost  effectiveness  with  which  coatings  can  be  applied.  These 
advantages  can,  however,  readily  be  lost  if  the  proper  chemistry  and 
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processing  conditions  are  not  cnosen.  For  example,  if  the  hydrolysis  and 
condensation  rate  of  one  of  the  components  in  a  multicomponent  system  is 
mucn  faster  than  the  others,  one  oxide  or  hydroxide  may  precipitate  out  of 
the  solution  (upon  the  addition  of  water)  before  the  other  has  begun  to 
react . 

Data  of  several  investigators  have  indicated  that  in  the  case  of 
tetraethylorthosilicate  (TEOS),  the  rate  of  disappearance  of  water  can 
be  described  with  simple  first  order  kinetics.  Assuming  this  to  be  the 
case,  the  absolute  reaction  rate  was  precisely  determined.  Rates  for  other 
alkoxides  could  theoretically  be  measured  and  compared,  and  calculations 
could  be  made  as  to  how  much  water  should  be  added  and  when  it  should  be 
added  to  maximize  the  number  of  bridges.  Here  and  are  metal 

cations  (Si,  Al,  etc.)  and  0  is  oxygen.  Pursuing  this  lead,  investigations 
were  undertaken  of  the  kinetic  behavior  of  TEOS  in  ethanol  and  water  using 
infrared  spectroscopy  to  measure  the  concentration  of  water  as  a  function 
of  time,  using  the  1640  cm  ^  band. 

All  the  infra-red  data  were  taken  on  a  Perkin  Elmer  model  297  using  a 
calcium  fluoride  cell  of  nominal  pathlength  0.01  mm.  The  absolute  path- 
length  was  not  calibrated  before  each  run,  so  that  absorbance  cannot  be 
converted  directly  to  concentration.  As  a  result,  the  data  cannot  be 
compared  from  run  to  run  except  where  the  cell  was  not  taken  apart  between. 
However,  each  experiment  may  be  analyzed  on  a  relative  concentration  basis. 
Fortunately,  in  determining  the  order  of  the  reaction,  absolute  concentration 
does  not  matter,  only  the  relative  rate  of  change  as  discussed  below. 

The  form  of  the  data  for  all  starting  compositions  appeared  similar 
when  the  log  of  the  relative  rate  was  plotted  against  the  log  of  the 
absorbance.  Figure  6  shows  a  representative  plot  for  the  1:1 : 4-starting 
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composition  (the  numbers  represent  the  molar  ratios  of  water  to  TEOS  to 
ethanol).  As  was  the  case  for  all  experiments  run,  the  acid  concentration 
of  the  entire  solution  was  10  “  normal  HC1 .  When  these  data  are  compared 
with  those  obtained  on  the  1:2:4  and  1:4:4  starting  compositions  (Figures 
7  and  8) ,  it  can  be  seen  that  the  initial  apparent  orders  of  the  reaction 
decrease  with  increasing  water  concentration  and  even  become  negative  in 
the  1:4:4  case.  All  three  systems  show  significant  deviation  from  first 
order  kinetics  throughout  most  of  the  course  of  the  reaction. 

In  an  effort  to  determine  if  starting  with  very  clean  components 
alters  the  apparent  order  of  the  kinetics,  several  experiments  on  the 
1:1:4  starting  composition  were  made.  As  illustrated  in  Figure  9a,  there 
was  little  difference  between  TEOS  used  as-received  compared  with  the  same 
material  which  had  been  freshly  distilled,  when  both  were  hydrolysed  with 
deionized  water  and  Acculute  HC1  was  used  as  the  acid  source.  There  was  no 
noticeable  difference  when  18  megaohm  water  was  used  for  hydrolysis  in 
conjunction  with  Acculute.  Finally,  silicon  tetrachloride  was  added  to 
ethanol  to  generate  the  HC1;  and  when  this  was  used  in  conjunction  with  18 
megaohm  water,  the  initial  apparent  order  fell  much  closer  to  unity 
(Figure  9b).  This  experiment  was  run  twice;  and  the  results  point  out  how 
much  variation  might  be  expected  in  such  data.  This  entire  series  of 
experiments  was  run  without  dismantling  the  cell,  so  that  the  results  can  be 
compared  quantitatively. 

Data  obtained  under  a  range  of  conditions  were  analyzed  in  terms  of 
the  principles  of  chemical  kinetics.  In  these  terms,  the  rate  of  disappearance 
of  one  of  the  species  in  a  reaction  can  be  related  to  the  concentrations  of 
the  reactants.  For  the  reaction  shown  as  Equation  (3),  the  rate  of 
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disappearance  of  A  can  be  expressed  as  Equation  (4) 

A  +  B  -  C  +  D  (3) 

-  ^  =  Kf[A]m[B]n  (4) 

Here  m  and  n  are  the  orders  of  the  reaction  with  respect  to  each  component. 
An  important  assumption  made  here  is  that  the  reaction  is  irreversible.  If 
this  were  not  the  case,  the  rate  of  production  of  A  using  the  back  reaction 
can  be  shown  as  Equation  (5)  : 

=  Kb[C]k[D]1  (5) 

The  total  rate  of  change  of  the  concentration  is  the  sum  of  these  two 
equations.  In  these  equations,  and  are  reaction  rate  constants  for 
the  forward  and  back  reactions,  respectively.  They  are  independent  of 
composition,  but  are  functions  of  temperature. 

In  the  standard  analysis,  the  concentration  of  A  is  plotted  against 
time,  and  a  polynomial  curve  is  fit  to  the  data.  The  derivative  is  taken, 
and  the  log  of  this  is  plotted  against  the  log  of  the  concentration.  In 
most  cases,  a  straight  line  results  whose  slope  is  equal  to  m,  the  order 
of  the  reaction  with  respect  to  A.  This  is  only  true,  however,  if  the 
concentration  of  B  is  relatively  constant.  If  this  is  not  the  case,  a  more 
complex  form  is  followed. 

Reaction  rate  data  are  usually  taken  with  all  components  except  one 
in  vast  excess  in  order  to  eliminate  the  aforementioned  problem.  In  our 
case,  however,  we  are  working  in  a  composition  range  where  the  initial 
concentrations  of  A  and  B  are  comparable.  In  the  case  where  the  initial 
concentrations  are  identical.  Equation  (4)  can  be  reduced  to  Equation  (6).: 


35. 


d  [  A] 
dc 


Kf  [A] 


nrt-n 


(6) 


If  the  reaction  is  actually  first  order  in  both  components,  it  will  appear 
second  order  with  respect  to  A  (due  to  B  varying  at  the  same  rate  as  A) . 

For  the  sake  of  simplicity,  the  orders  of  reaction  with  respect  to  each 
component  will  be  assumed  to  be  unity  from  this  point  on,  as  is  assumed  to 
be  the  case  for  TEOS. 

Upon  separation  of  the  variables  and  integrating,  the  concentration  as 
a  function  of  time  is  given  as  Equation  (7): 


[A]  = 


[A], 


Kct [ A1  -  1 
f  r 


(7) 


For  an  excess  of  B,  the  derivation  is  slightly  more  difficult,  and  has  been 
outlined  below  for  the  initial  concentrations  [B]  *  4 [A] . 

Equation  (3)  shows  that  A  and  B  combine  in  equimolar  quantities,  so 
that  for  the  initial  conditions  of  [B]  =  4 [A],  there  is  a  threefold  excess 
of  B.  If  x  is  the  concentration  of  A  as  a  function  of  time,  then  the  con¬ 
centration  of  B  equals  x  +  3.  Substituting  this  expression  into  Equation 
(4),  Equation  (8)  is  generated. 


dx 

dt 


*  K^x(3+x) 


(8) 


Separating  the  variables,  one  obtains  Equation  (9): 


dx 

x(3+x) 


-  Kfdt 


(9) 


Setting  up  the  definite  integrals  of  x  going  from  [Aq]  to  [A]  and  time 
going  from  0  to  t,  one  obtains  Equation  (10): 


in  <^J)  -  ln(i)l 


t  - 


(10) 


36. 

Following  the  same  procedure  for  the  general  initial  case  [B]  =  [A]+y,  one 
obtains  the  general  form: 

[A] (y+[ A  ] ) 

t  VK  ^-n  (11) 

-VKf  tAo](y+[A]) 

In  order  to  examine  the  apparent  order  of  the  reaction  as  evidenced  by 
the  disappearance  of  water,  it  is  necessary  to  plot  the  log  of  the  rate  of 
disappearance  against  the  log  of  the  concentration  and  measure  the  slope. 

The  method  described  below  is  roundabout,  but  parallels  the  analysis  of 
the  computer  simulations  of  hydrolysis  which  we  have  carried  out.  Using 
Equation  (11),  a  table  of  water  concentration  vs.  time  was  generated.  These 
data  were  then  normalized  by  scaling  the  initial  concentration  to  1  and 
choosing  values  for  so  that  all  curves  passed  through  another  common 
point.  This  procedure  makes  the  A  to  B  ratio  the  only  variable. 

At  this  point,  a  plot  of  the  normalized  concentration  vs.  normalized 
time  was  generated  (Figure  10) .  The  derivatives  of  the  curves  were  cal¬ 
culated  and  ln(-  -^^-)  was  plotted  against  ln[A]  (Figure  11).  Several  features 
of  this  graph  should  be  noted.  An  initial  ratio  of  B  to  A  of  16  is  needed 
for  the  initial  part  of  the  reaction  to  appear  to  be  first  order  in  A.  As 
this  ratio  decreases,  the  apparent  order  increases.  However,  the  initial 
apparent  order  increases  faster  than  the  apparent  order  during  the  latter 
stages  of  the  reaction.  This  results  in  a  curve  that  is  initially  steep, 
and  approaches  a  slope  of  unity  as  the  reaction  proceeds.  However,  when 
A  and  B  are  present  in  stoichiometric  amounts,  the  slope  is  constant  at  2. 

This  shows  that  the  apparent  order  for  even  the  initial  part  of  the  reaction 
can  vary  from  one  to  two,  depending  on  the  initial  concentration  rate  of  A  to  B. 
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Up  to  now,  we  have  discussed  a  system  which  undergoes  only  one  possible 
reaction  which  is  first  order  and  irreversible.  In  the  case  of  TEOS ,  these 
assumptions  are  incorrect.  Water  is  a  reactant  in  all  four  of  the  possible 
hydrolysis  reactions  and  half  of  the  many  condensation  reactions,  so  even 
if  we  maintain  the  irreversibility  criterion,  the  behavior  of  the  water  con¬ 
centration  as  a  function  of  time  may  be  quite  complex. 

Consider  only  two  possible  reactions  which  use  a  common  reactant,  as 
indicated  in  Equations  (12)  and  (13): 


A  +  B  -*■  C  +  D 

(12) 

A  +  C  E  +  D 

(13) 

These  two  reactions  correspond  to  the  first  two  hydrolysis  reactions 
of  TEOS,  where  A  is  water,  B  is  TEOS,  C  is  triethoxysilanol ,  D  is  ethanol, 
and  E  is  diethoxysilanediol .  If  we  were  able  to  inhibit  one  of  these 
reactions  while  allowing  the  other  to  proceed,  we  would  find  first  order 
kinetics  for  the  disappearance  of  water  when  we  start  with  a  vast  excess  of 
silane.  If  we  were  to  assume  that  the  rate  of  the  second  reaction  is  higher 
than  the  first  (as  is  predicted  for  TEOS)  while  still  occurring  independently, 
we  would  expect  results  similar  to  Figure  12. 

For  the  case  of  TEOS,  these  two  reactions  are  not  separable.  However, 
when  the  reactions  are  initiated,  since  reactants  for  only  the  first  reaction 
are  present,  the  instantaneous  rate  of  the  reaction  would  lie  on  the  lower 
of  the  two  lines  in  Figure  7.  As  soon  as  the  first  product  appears,  the  rate 
of  water  consumption  may  deviate  from  this  line.  For  the  case  where  B  is 
present  in  vast  excess,  if  the  rate  of  the  second  reaction  is  infinitely 
faster  than  the  first,  we  would  expect  to  see  an  overall  rate  of  water 
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consumption  double  that  of  the  first  reaction.  This  is  because  as  soon  as 
a  molecule  of  B  reacts,  a  molecule  of  C  will  be  produced  whicn  will  instan¬ 
taneously  consume  another  water  molecule.  This  will  be  reflected  by  a  line 
that  still  has  a  slope  of  unity,  but  is  displaced  upward  by  a  factor  of 
ln(2).  In  the  other  extreme,  we  have  already  said  that  if  C  does  not  react, 
the  overall  rate  will  follow  the  lower  line. 

It  has  been  mentioned  earlier  that  at  least  a  sixteen-fold  excess  of 
B  is  needed  for  the  slope  of  the  ln(-rate)  vs.  In (concentration)  to  be  unity. 
If  the  previous  analysis  is  carried  out  in  a  regime  where  the  starting  B 
to  A  ratio  is  only  8,  the  results  change.  As  we  have  seen  before,  if  C 
does  not  react,  the  slope  of  the  initial  curve  will  be  slightly  greater  than 
unity.  If  the  reaction  rates  were  the  same,  we  would  have  a  situation 
equivalent  to  having  a  sixteen-fold  excess  of  B  where  C  does  not  react.  And, 
as  before,  if  C  reacts  instantaneously  we  would  expect  a  curve  parallel 
to  the  original  one,  displaced  higher  by  a  factor  of  ln(2).  For  the  case 
where  C  reacts  faster  than  B  but  is  not  instantaneous,  we  would  expect  the 
resulting  overall  rate  to  fall  between  the  two  latter  cases  (Figure  13). 

These  predictions  are  born  out  in  the  computer  simulation.  Figure  14  shows 
the  effect  of  increasing  the  rate  of  the  second  reaction  from  one  fifth  to 
ten  times  the  rate  of  the  first. 

In  addition  to  the  hydrolysis  reactions,  re-esterification  and  con¬ 
densation  also  involve  water,  here  as  a  product.  Re-esterification  is  by 
definition  the  reverse  of  the  hydrolysis  process  and  has  been  shown  to  occur 
at  significant  rates  for  silanes  in  studies  using  Karl-Fischer  reagent.  The 
result  of  reversibility  being  included  in  the  kinetic  treatment  is  to 
increase  gradually  the  apparent  order  of  the  overall  reaction  as  it  proceeds. 
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Figure  15  shows  Che  result  of  a  computer  simulation  of  re-esterification 
rates  between  0  and  10  percent  of  tne  forward  rate  of  one  reaction.  As  seen 
before,  when  the  hydrolysis  reaction  is  irreversible  and  the  silane  species 
is  present  in  great  excess,  the  resulting  apparent  order  is  unity.  However, 
as  the  rate  of  the  reverse  reaction  becomes  more  appreciable,  the  apparent 
order  increases  greatly  as  the  water  concentration  drops.  Finally,  as 
equilibrium  is  approached  the  overall  reaction  rate  goes  to  zero,  causing 
the  ln(-rate)  vs.  ln(conc.)  curve  to  tail  off  towards  negative  infinity.  As 
expected,  as  the  rate  of  the  re-esterification  process  increases,  the 
equilibrium  concentration  also  increases.  The  usual  approach  for  correcting 
this  effect  is  to  reassign  zero  as  the  equilibrium  concentration  of  the 
monitored  species.  This  serves  to  make  the  corrected  concentration  go  to 
zero  at  the  same  rate  as  the  reaction  rate,  preserving  the  first  order 
relation. 

Condensation  is  the  process  by  which  small  species  combine  to  form  larger 
ones.  It  can  be  represented  by  two  reactions,  one  involving  water  as  a 
product  (Equation  15)  and  the  other  yielding  alcohol  (Equation  16). 


sSi-OH  +  HO-Si=  -»• 

sSi-0-Si= 

+  H20 

(15) 

=Si-0H  +  RO-Si=  -*• 

*Si-0-Sis 

+  ROH 

(16) 

Since  the  latter  does  not  involve  water,  we  will  ignore  it  in  this  treatment. 
As  before,  water  produced  via  reaction  (15)  cannot  be  distinguished  from 
water  which  was  the  product  of  re-esterification.  However,  in  this  case, 
the  silane  product  cannot  be  fed  back  into  the  hydrolysis  reaction.  An 
extreme  example  of  this  would  be  the  hydrolysis  and  condensation  of 
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Figure  15 


tetramethylorthosilicate  (TMOS)  with  4  moles  of  water  per  mole  of  TMOS. 

Since  hydrolysis  of  TMOS  is  very  fast,  all  of  the  reactants  would  initially 
combine  to  form  silicic  acid,  which  would  then  condense.  All  four  moles  of 
water  would  disappear  followed  by  the  reappearance  of  two.  This  will  cause 
the  rate  of  disappearance  to  become  negative,  and  the  natural  log  will 
become  undefined.  In  TEOS,  hydrolysis  and  condensation  occur  simultaneously. 

An  added  complication  is  that  the  condensation  rate  will  be  controlled  by 
the  mobility  of  oligomeric  species  during  the  later  stages  of  condensation. 

These  processes  cannot  be  modelled  with  simple  kinetics. 

The  form  of  the  infra-red  data  obtained  in  the  present  work  cannot 
be  regarded  as  surprising  in  light  of  the  many  simultaneous  processes 
occurring.  Initial  apparent  orders  less  than  one  can  be  explained  by  subsequent 
steps  of  the  hydrolysis  being  faster  than  the  first.  As  the  reaction  proceeds, 
re-esterification  and  condensation  become  appreciable,  and  the  apparent  order 
will  steadily  increase.  These  explanations  will  account  for  all  of  the  trends 
seen. 

There  are  other  possible  explanations,  however.  If  we  assume  that  the 
apparent  order  is  in  fact  unity,  a  lesser  value  may  be  observed  for  the 
following  reasons:  Firstly,  the  temperature  of  the  fluid  inside  the  IR  cell 
may  increase  during  the  reaction  due  to  the  exothermic  nature  of  hydrolysis, 
or  heating  from  the  incident  infra-red  radiation.  Secondly,  the  activity 
of  the  hydrogen  ions  may  increase  over  the  course  of  the  reaction  as  the 
water  concentration  decreases  and  in  some  cases  goes  to  zero.  Thirdly,  the 
partially  hydrolysed  silane  may  exert  a  self-catalytic  effect,  increasing 
the  rate  of  the  reaction  as  it  proceeds.  The  first  possibility  has  seemingly 
been  ruled  out  by  workers  using  a  water-cooled  temperature-controlled  IR  cell. 
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and  by  ourselves  by  comparing  the  reaction  rates  with  varying  incident 
infra-red  intensities.  We  have  determined  that  for  our  simple  non-temperature 
controlled  calcium  fluoride  cell  with  a  0.5  mm  pathlength  (50  times  the 
experimental  conditions),  a  4°K  temperature  rise  will  occur  due  to  beam 
heating  of  ethanol.  This  result  must  be  taken  with  a  grain  of  salt  because 
the  volume  to  surface  ratio  has  been  increased  by  the  same  ratio,  decreasing 
the  rate  at  which  the  heat  can  be  conducted  away. 

The  possibility  involving  the  activity  of  hydrogen  ions  also  seems 
unlikely.  Our  data  indicate  that  the  rate  of  hydrolysis  decreases  if  the 
concentration  of  hydrogen  ions  remains  constant  but  more  water  is  added 
initially.  This  is  consistent  with  the  pKa  of  protonated  water  being 
higher  than  that  of  protonated  ethanol;  and  as  the  water  disappears  the 
hydrogen  ions  are  available  to  protonate  the  other  solvent  species,  i.e.  , 
ethanol.  Based  on  this,  the  effect  is  predicted  to  be  largest  in  systems 
where  the  water  just  disappears  completely.  This  is  contrary  to  what  we 
observed,  where  the  effect  is  most  acute  for  the  1:4:4  starting  composition. 

The  third  possible  explanation,  involving  an  autocatalytic  effect  of 
the  partially  hydrolysed  silane  has  no  supporting  data,  but  cannot  be 
ignored  as  a  possibility. 

In  summary,  we  have  measured  the  concentration  of  H^O  as  a  function  of 
time  during  the  hydrolysis  of  TEOS  using  IR  spectroscopy.  The  apparent  order 
of  the  reaction,  as  calculated  from  the  water  concentration,  changed  as  the 
reaction  proceeded.  It  is  our  opinion  that  this  is  due  in  large  part  to 
following  only  one  species,  water,  which  is  involved  in  several  hydrolysis 
and  condensation  reactions.  As  each  of  these  reactions  likely  takes  place 
at  different  rates,  the  resulting  convolution  makes  the  data  complex. 
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In  order  to  overcome  this,  there  are  two  choices.  The  first  is  to 

measure  only  the  water  concentration  but  to  monitor  only  the  start  of  the 

reactions  for  very  understoichiometric  amounts  of  water,  far  below  the 

conditions  discussed  here  and  of  interest  in  the  development  of  coatings. 

Although  this  may  result  in  a  quantitative  reaction  rate  and  order,  there 

is  some  doubt  to  how  applicable  these  numbers  would  be  to  systems  with 

more  water.  TEOS  does  not  obey  simple  kinetics,  so  that  there  would  be 

much  difficulty  in  extrapolating  the  data  to  higher  water  systems. 

Alternatively,  one  could  follow  the  concentration  of  all  species  as  functions 
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of  time.  This  may  be  possible  with  tools  such  as  Si  NMR  if  done  along 
the  same  lines  as  condensation  studies.  Although  the  ensuing  analysis 
would  be  difficult,  it  does  not  seem  impossible. 

Infra-red  spectroscopy  does  not  lend  itself  well  to  this  task. 

Absorption  bands  are  generally  quite  broad  and  tend  to  overlap,  leading 
to  difficulties  in  determining  how  much  of  any  species  is  actually  present. 

IR  techniques  do  make  the  determination  of  water  concentration  quite  easy; 
and  the  point  at  which  the  free  water  in  solution  disappears  can  be 
determined  without  difficulty.  Such  data  are  of  great  utility  when 
partially  hydrolysed  silanes  are  added  to  other  faster  reacting  species. 

In  the  past,  one  has  relied  upon  experience  to  decide  how  long  after  the 
silane  has  started  hydrolysing  to  add  other  species.  IR  analysis  makes 
this  determination  quite  easy  and  allows  the  maximization  of  M-O-Si 


bridges. 
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B.  Crystallization  vs.  Viscous  Sintering 

Using  soi-gel  techniques,  it  is  possible  to  produce  fine  amorphous 
powders  or  dried  gels  having  distributions  of  fine  pores.  To  obtain  dense, 
glassy  bodies  or  coatings  from  such  materials,  it  is  useful  to  define 
regions  of  time  and  temperature  (or  heating  rates)  where  the  kinetics  of 
viscous  sintering  are  faster  than  those  of  crystallization — assuming  that 
such  regions  exist  for  a  given  material. 

In  describing  viscous  sintering,  we  have  used  the  cylindrical  pore 
model  of  Scherer.  This  model  is  very  attractive  for  bodies  of  the  type 
produced  by  sol-gel  processes.  It  employs  the  basic  assumptions  of  the 
Frenkel  model  for  the  viscous  sintering  of  spheres;  and  considering  these 
assumptions,  it  is  remarkable  that  the  model  predicts  experimental  data  as 
well  as  it  does  (which  is  quite  well). 

The  model  predicts  that  the  time  required  for  appreciable  sintering  is 
proportional  to  the  viscosity  and  inversely  proportional  to  the  surface 
energy.  Among  microstructural  features,  the  most  important  feature  is  the 
pore  radius.  The  initial  particle  size  and  packing  density  are  less 
significant  (particles  of  difference  sizes  can  be  packed  differently  to  give 
the  same  pore  size,  and  exhibit  similar  times  for  sintering).  In  the 
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present  paper,  we  have  taken  pore  sizes  of  50  A  and  5000  A  to  represent 
the  effective  limits  obtained  from  sol-gel  processing. 

The  process  of  crystallization  has  been  described  using  the  approach 
developed  previously  in  our  laboratory.  Under  isothermal  conditions,  the 
times  at  various  temperatures  required  to  reach  a  just-detectable  degree 
of  crystallinity  are  represented  by  time-temperature-transformation  (TTT) 
curves;  while  crystallization  under  conditions  of  constant  heating  rate  is 
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described  using  the  analysis  of  crystallization  statistics.  In  the  present 
paper,  viscous  sintering  and  crystallization  are  viewed  as  competitive 
processes  and  the  analysis  of  this  competition  is  applied  to  a  number  of 
representative  glass-forming  liquids. 

Considering  isothermal  conditions,  Figure  17  indicates  the  calculated 
behavior  for  a  representative  organic  glass  former,  o-terphenyl.  The 
regions  of  time  at  any  temperature  between  the  appropriate  sintering  curve 
and  the  crystallization  curve  for  a  just-detectable  volume  fraction 
crystallized  at  10  ^  indicate  the  conditions  where  densif ication  can  be 
achieved  without  detectable  crystallization.  In  some  cases,  the  regions 
of  time  may  be  inaccessibly  short;  but  in  the  case  of  o-terphenyl,  there 
is  a  range  of  conditions  under  which  densif ication  can  be  obtained  without 
crystallization. 

The  magnitude  of  the  viscosity  at  any  temperature  affects  both  the 
sintering  rate  and  the  crystallization  rate.  To  take  account  of  the 
dependence  of  the  viscosity,  results  such  as  those  shown  in  Figure  17 
may  conveniently  be  plotted  as  a  function  of  t/n,  where  n  is  the  viscosity. 
Such  a  plot  is  shown  in  Figure  18  for  anorthite. 

For  crystallization  under  conditions  of  constant  heating  rate,  the 
competition  between  crystallization  and  viscous  sintering  is  conveniently 
represented  in  terms  of  a  heating  rate  vs.  temperature  relation,  such  as 
that  shown  in  Figure  19  for  anorthite.  As  shown  there,  there  exist 
windows  in  temperature,  over  a  range  of  heating  rates,  in  which  substantial 
sintering  can  take  place  without  the  occurrence  of  detectable  crystllization. 
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All  of  the  above  discussion  has  used  bulk  material  properties  and 
measured  or  estimated  nucleation  barriers,  together  with  the  assumption 
of  steady  state  nucleation  rates.  The  effect  of  different  nucleation 
barriers  on  the  competition  between  sintering  and  crystallization  is 
shown  in  Figure  20  for  anorthite-like  materials  having  various  nucleation 
barriers  (represented  as  BkT).  As  seen  there,  with  decreasing  barrier  to 
crystal  nucleation,  the  processing  window  in  which  sintering  can  be 
achieved  without  detectable  crystallization  is  expected  to  decrease.  For 
sufficiently  small  nucleation  barriers,  such  as  BkT  =  40  for  a  material 
with  the  properties  of  anorthite,  crystallization  will  take  place  before 
sintering  at  all  heating  rates  evaluated.  The  measured  nucleation  barrier 
of  anorthite  is  about  82  kT. 

The  effect  of  changes  in  the  degree  of  crystallinity  is  shown  in 
Figure  21  for  Na20*2Si02.  Larger  processing  windows  are  suggested  for 
conditions  where  larger  volume  fractions  crystallized  are  acceptable.  Also 
plotted  in  Figure  21  are  the  experimental  results  of  Hench  et  al.,  who 
observed  sintering  without  crystallization  in  one  hour  treatments  at 
753  C  -  793  C,  and  crystallization  at  temperatures  of  793  C  -  813  C. 

Their  condition  for  sintering  corresponds  to  a  degree  of  densification 
beyond  that  contemplated  in  the  present  calculations,  and  hence  the 
agreement  between  predictions  and  experiment  seems  quite  reasonable. 

The  assumption  of  steady  state  nucleation  rates  can  be  relaxed  by 
considering  the  incubation  time  during  which  steady  state  distributions  of 
subcritical  embryos  are  established.  The  effect  of  such  transient 
nucleation  is  to  shift  the  crystallization  curves  to  higher  temperatures 
and  thereby  enlarge  the  processing  window.  Of  the  materials  considered 
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in  the  present  study,  the  effect  is  most  pronounced  for  L^O-ZSiO  ,  for 
which  the  results  shown  in  Figure  22  were  obtained. 

It  is  well  recognized  that  sol-gel  derived  materials  are  generally 
characterized  by  properties,  such  as  liquid  viscosities,  which  differ 
from  those  of  melted  glasses.  In  applying  the  present  analysis  to 
SiC^,  use  has  been  made  of  the  data  of  Sacks  et  al.  for  the  viscosity 
of  the  "wet"  SiO^  employed  in  his  studies.  The  viscosities  reported  there 
were  appreciably  smaller  than  those  of  conventional  fused  silica,  and  the 
lower  viscosities  have  an  appreciable  effect  on  the  expected  rates  of 
sintering  and  crystallization.  These  are  shown  in  Figure  23.  Also  shown 
in  the  figure  are  the  other  results  of  Sacks,  who  observed  sintering  to  an 
amorphous  state  in  24  hours  at  1274  C,  but  crystallization  in  48  hours 
at  1273  C  or  3  hours  at  1323  C.  The  pore  size  in  Sacks'  experiments  were 
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a  few  thousand  A.  From  the  results  shown  in  Figure  23,  remarkable  agreement 
is  seen  between  the  predicted  behavior  and  experimental  results. 

The  effects  of  other  factors  on  the  competition  between  crystallization 
and  viscous  sintering  are  presently  being  explored.  These  include  the 
effects  of  heterogeneous  nucleation  as  well  as  the  effects  of  time-dependent 
properties.  Once  the  results  of  these  calculations  are  in  hand,  it  should 
be  possible  to  specify  for  a  given  material  the  conditions  of  time  and 
temperature  in  which  densification  of  coatings  can  be  achieved  without 
undesired  crystallization  (or  alternatively,  the  conditions  under  which 
crystallization  to  desired  assemblages  of  phases  can  be  achieved) .  This  will 
greatly  facilitate  the  experimental  development  of  coatings  for  glass  having 
the  desired  high  density  (freedom  from  porosity).  It  is  expected  that  the 
theoretical  work  will  be  completed  by  the  time  that  the  optimization  phase 
of  the  respective  coating  developments  is  initiated. 
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